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ABSTRACT 



We pre sent an updated ana l ysis of the M31 pixel lensing candidate event 0AB-N2 previously re- 
ported in ^ Calchi Novati et alj (l2009a). Here we take advantage of new data both astrometric and 
photometric. Astrometry: using archival 4ni-KPN0 and HST/WFPC2 data we perform a detailed 
■ analysis on the event source whose result, although not fully conclusive on the source magnitude 

determination, is confirmed by the following light curve photometry analysis. Photometry: first, 
Q , unpublished WeCAPP data allows us to confirm 0AB-N2, previously reported only as a viable can- 

$—( ' didate, as a well constrained pixel lensing event. Second, this photometry enables a detailed analysis 

' in the event parameter space including the effects due to finite source size. The combined results of 

, these analyses allow us to put a strong lower limit on the lens proper motion. This outcome favors 

the MACHO lensing hypothesis over self lensing for this individual event and points the way toward 
distinguishing between the MACHO and self-lensing hypotheses from larger data sets. 
^ ■ Subject headings: dark matter — gravitational lensing — galaxies: halos — galaxies: individual (M31, 

V/-^ ■ NGC 224) — Galaxy: halo 

00 

^ , 1. INTRODUCTION 

\ Gravitational microlensing is an established tool of research. The original motivation ()Pacz vnskill 1 986[ ) has been the 
. search for dark matter in the form of massive compact halo obje cts (MACHOs) (for a recent discussion from a larger 
^—^ ■ point of view on dark matter and gravitational lensing we refer to lMassev et al.ir2010[ ). The first lines of sight explored 
to this purpose have been those toward the Magellanic Clouds (LMC and SMC). M icrolensing campaigns toward the 
Galactic center have then been used to constrain Galactic st ructure (iMonied 120101 and references therein) and, more 

recently, for the detection of extra-solar planets (jCaudi 20 13; lDominildl2010l and references therein) . 

rS Results o n MACHOs toward t he LMC and the SMC have b een reported by the MACHO ( Alcock et all l200(/). 
the EROS (jTisserand et al.ll2007l ) and the OGLE collaboration ()Wvrzvkowski et al.l 12009 1. There is agreement that 
MACHOs are excluded as a major component of dark matter for a wide range of possible mass values (down to 
about 10~^ M(7)) . There remain s however a range, (0.1 — 1) M0, in which the results of the MACHO collaboration, 
as confirmed bv iBennettI (|2005| ). indicate that MACHOs might compose about a / = 20% fraction of the Milky 
Way (hereafter MW) halo mass at least up to the distance of the LMC. This outcome is challenged by the results 
of the EROS and the OGLE collaborations. A possible contamination of the microlensing signal due to MACHOs, 
once the background of intrinsic variable stars that can mimic microlensing is excluded, is self lensing, w herein both 
the source and lens belong to some luminous population^ . (This possib ility has be en first conside red in ISahul 119941 : 
[Gould 19 9^ and further s tudied by several authors, e. g. iGvuk et"aL 2000; Ma ncini et al.l 120041: we refer also to 
[Calchi No vati et al.ll2009bl for a specific analysis of the OGLE-II results, to Calc hi Novati et al.ll200 6 for a study of the 
possible contribution of the LMC own dark matter halo to the observed MACHO events and to|Moniez 2010 for an 
updated review on this issue). Indeed, the fact that the only remaining allowed mass range for MACHOs, from the 
LMC analyses, corresponds to that of the stars that might also act as lenses may be indicative of some observational 
bias (we recall that the characteristics of the microlensing events, in particular their duration, depend on the lens 

* Royal Society University Research Fellow 
t http:/ /plan. physics. unisa.it/index. htm 

1 Hereafter we will refer to "self lensing" , broadly speaking, for any microlensing event with a stellar lens, not necessarily belonging to 
the same stellar population of the source. 
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mass). On the other hand, if the reported events are reahy to be attributed to MACHOs, a fortiori for such a sizeable 
fraction as that imphed by the MACHO cohaboration results, this might have some deeper astrophysical meaning. 

Beyond the Magell anic Clouds, t he next sui table ta rget for a mi crolensing search is our nearby and similar galaxy 
of Andromeda, M31 (jCrottsI [19921 iBaillon et~ al. 1993; lJetzeij[T994l) . It allows one to explore a different line of sight 
through the MW halo; we can fully map MBl's own dark matter halo (which is not possible for the MW, this 
being perhaps the most severe limitation of LMC studies); finally, the inclination of M31 is expected to give rise 
to a characteristic signature in the spatial distribution of M31 halo events, so facilitating their identification against 
the contamination of self-lensing events. Because of the distance of M31, the sources for microlensing eve nts are 
no lo nger resolved objects, at least for ground-based telescopes, so that we enter the pixel lensing regime (IGouldl 
19961). Several observational campaigns have already been undertaken along this line of sight (for a review we refer to 
Calchi Novati|[2009) . The resuhs obtained up to now on the MACHO issue are controversial. The POINT- AGAPE 
collaboration claimed evid ence of a MACHO signal in the same mass range indicated by the MACHO LMC analysis 
(jCalchi Novati et al.ll2005l ). This outcome, however, has been challenged by the MEGA colla boration who concluded 
that the detected signal is compatible with the expected self-lensing rate ()de Jong et al.|[2006[) . Finally, the WeCAPP 
collaboration is expected to present the final analysis of their 11-year campaign^. 

The issue of self lensing for M31 microlensing is even more difficult than for the LMC for two principal reasons. 
First, the expected self-lensing rate is quite large with respect to MACHO lensing, at least in the inner part of M31 
(in which both signals have the larger expected rate because of the huge number of available sources). The exact ratio 
depends not only on the unknown halo fraction in the form of MACHOs and on the MACHO mass but also on the not 
fully understood contribution of M31 stellar lens ing. (In f act, a possib le approa ch to deal with this problem is that 
followed in the already cited analyses of Calchi Novati et al. 2005 and de Jong et al. 2006. where statistical arguments 
were used for a full set of events. Interestingly, a main point of disagreeme nt between these two analy ses is the role 
played by self lensing with respect to MACHO lensing, as discussed also in llngrosso et al.ll2006L 120071 ). Second, the 
lack of knowledge of the source flux, specific for the "pixel" lensing regime, adds a further degeneracy in the event 
parameter space that makes the disentangling of the two signals more difficult. In particular, this degeneracy implies 
that, rather than determining the Einstein time, ^e, the impact parameter, mq, and the source flux , from the ligh t 
curve one can usually reliably estimate only the full- width- half- maximum duration, tpwHM = /(wq) ()Gondololll999t l. 
and the flux devi ation at maximum amplifica tion. At/) = (jj* Amax('Uo), usually expressed in term of magnitude, e. g. 
Ai? (|Gouldlll99fil : iWozniak fc Paczvnskil 119971 ). 

This problem of interpretation of the lens nature, together with the overall small number of reliable microlens- 
ing candidate events reported toward M31 so far, make the issue of a correct astrophysical understanding of single 
events that may contain precious additional infor mation particularly re levant. For M31 lensing this issue has been 
first emphasized by the WeCAPP collaboration in iRiffeser et al.l (|2006l ). The WeCAPP collaboration, furthermore, 
presented an extremely detai l ed analysis of the mic rolensing event P A-S3/GL1 (first rep orted by POINT-AGAPE, 
iPaulin-Hen riksson et al.|[2003t iBelokurov et al.ll2005l . and subsequently, IRiffeser et al.|[2003l detected in their own data 
set by the WeCAPP collaboration itself). Also thanks to the excellent sampling al ong the bump assu red by the two 
independent data sets, and on the basis of a study of the source finite size effect (iWitt fc Maolll994D . IRiffeser et al.l 
(|2008i) showed that this extremely bright microlensing event is more likely to be attributed to MACHO lensing than 
to self lensing. Interestingly, this would hold even though the event is located very near to the M31 center, d ^ 4', 
where the expected self-lensing rate is larger. 

In fact, the first detailed analysis of a M31 pixel lensing event aimed at t he characterization of the lens nature 
has been presented by the POINT- AG APE cohaboration for the event PA-Nl (jAuriere et al.ll2001f ). With the source 
magnitude and color fixed thanks to an analysis of archival Hubble Space Telescope (HST) data, and based on a specific 
Monte Carlo simulation of the experiment, where the observed characteristics of the evei it had been us e d to d elimit the 
possible parameter space (using in particular a lower limit for the lens proper motion) , lAuriere et al.l ()2001h conclude 
that the PA-Nl lens is more likely a MACHO if the mass halo fraction in the form of compact halo object is above 
20%. The results of thi s analysis have been however challen ged by the MEGA colla boration (who also reported this 
event as MEGA-ML16. Ide Jong et al]l2006l ). in particular bv lCseresnies et all (I2005D who argued, on the basis of new 
HST observations, against the source identification reported in lAuriere et al.l ()200lD . 

Following the strong interest in the characterization of single events, in the present work we will present an updated 
analysis o f the microlensing candidate event 0AB-N2 first reported by the PLAN (Pixel Lensing ANdromeda) collab- 
oration in lCalchi Novati et al.l ([2009a). First, we present a new analysis related to the issue of the source identification 
and characterization using KPNO and HST archival images. Sect. 12.11 The driving motivation for a new analysis 
comes from the significa ntly improved samp ling along the flux deviation due to the access to unpublished data of the 
WeCAPP cohaboration (jRiffeser et al.ll2001[ ). As we show, these new data allow us to firmly establish the microlensing 
nature of this flux variation and to carry out an improved and more detailed analysis of the lensing parameter space 
with the inclus fon of the fi nite source size effect. This allows us to put some constraints on the relative lens-source 
proper motion (|Gouldlll99l ). Sect.O In Sect. [3] we discuss the outcomes of the analysis within the framework of the 
Monte Carlo simulation of our M31 pixel lensing experiment. The limits on the lens proper motion are here relevant 
as they can be used as a tool to distinguish bewteen the self lensing and MACHO lensing nature of this microlensing 
event. 

^ A. Riffeser and S. Seitz, private communication. 
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2. ANALYSIS 
2.1. Source identification and characterization 

In iCalchi Novati et al.l (|2009aD we have tent atively identifi e d the 0AB-N2 candidate event source with a rather 
bright star, R = 20.9 R — I = 1.2, the nearest fMassev et al.l (|2006f ) catalog object to our estimated event position. 
Also motivated by the updated light curve analysis that suggests a fainter source star (Sect. 12. 2p we have performed 
a more detailed analysis. As a result we find this catalog star not to be the event source. However, the available data 
do not allow us to conclusively identify and characterize the candidate event source. 

More specifically, we have proceeded as follows. We have downloaded the images used by iMassev et al.l ()2006t ) for 
their analysis, collected at the 4m Mayall-KPNO telescope"^. In particular we have used the / band image M31F5I- 
l.fits with exposure time of 150 sec. We have selected a sample of ^ 200 stars brighter than Rc = 20 identified 
both on M31F5I-1 and our OAB image to carry out a pixel-to-pixel relative astrometry calibration^, for which we get 
to rms ~ 0.2". Because at maximum amplification the event is a resolved object, the accuracy of the event position 
determination is better, below 0.1". Taking therefore a = 0.2" we find that the catalog star previously identified as 
the event source is in fact about 4.3a away from the event position (and this is also confirmed to be the star reported 
in the catalog nearest to the eve nt position), wherea s a fainter object, not present in the catalog but still clearly visible 
(as checked also on a few other iMassev et aLl '2006 images), is located at only 1.3(7 (hereafter we refer to the catalog 
object as SI and to the fainter one as S2). This result is confirmed by an independent a nalysis in which bo th the 
KPNO and the OAB images are first "astrometrised" with respect to the USNO-Bl catalog (jMonet et al.ll2003[) . This 
analysis, however, has worse precision since the nearest catalog objects that avoids the inner M31 central region are 
located at about 3' — 4' from the event position. 

As a second step we perform PSF photometry analysis of the KPNO /-band image. For S2 we estimate Ic — 
20.5 ± 0.2. As expected, this is fainter than SI, by about ~ 1 mag, but the statistical error turns out to be extremely 
large, more than twice greater than the average error we find for stars of similar brightness in this image. Furthermore, 
we find also the output fit sharpness parameter for S2 to be a clear outlier. Both these elements indicate that S2 is 
more likely to be a blend of two or more stars. 

To investigate further this issue we have looked for archival HST images of this field. We find 2 F606W and 3 F814W 
WFPC2/WF3 images (pixel scale 0.1") available (GO 5971, P.I. Griffiths), out of which, however, only a 1000 sec 
exposure time F814W image is reported not to have quality problems, which we then use for our analysis. The cross- 
match of the KPNO image with the HST one is carried out by the identification of 12 stars on both images and we 
end up with a relative astrometry precision of rms ^ 0.1". Corresponding to the KPNO S2 position we find two stars 
(as identified also in the HST photometry analysis) lying within a distance of about 0.3" (corresponding roughly to 
one KPNO pixel). This confirms our previous conclusion on S2 being a blend. Both of these stars lie within 3ct of 
the estimated event position, with the brighter one lying nearer (about la). Finally, we have performed a photometry 
analysis on the HST image. In particular w e have made use of both the HST Stetson^ PSF and of the HSTphot 
WFPC2 photometry package ()Dolphinll2000[ ). for which we have considered three options for the background fit and 
PSF evaluation. In accord with this latter analysis we find, for the brighter of the two sources lying at the S2 position, 
-^instr ~ 20.6 — 21.1^, depending on the photometry option, with a formal statistical error a ~ 0.1 mag, and a magnitude 
difference between the two sources of ^ 0.2 mag. For the first analysis we find a somewhat brighter magnitude value 
for the brighter object, although with larger error, a ~ 0.2 — 0.3 mag, and larger magnitude difference, ~ 0.8 mag. 
Taking the HSTphot analysis as the fiducial one, we get to a final estimate of /instr = 20.8 ± 0.3 for the brighter star. 
The lack of color information, the large M31 background level and crowdedness at the event position which makes the 
photometry less reliable, and the presence of a second, fainter, star as a viable source do not allow us to consider these 
results on the HST photometry as conclusive. 

Overall, this new astrometry and jmaf/e-based photometry analysis, although with the already stressed caveats, 
allows us to definitively affirm that the 0AB-N2 source is a fainter star than previously concluded. This result turns 
out to be in agreement with the outcome of the new photometry light cwrwe-based analysis that we will detail in the 
following Section. 

2.2. Light curve analysis 

The microlensing flux variation 0AB-N2 has been first reported in ICalchi Novati et al.l (|2009al ) in the frame- 
work of a long-term pixe l lensing monitoring of M31 by the PLAN collaboration at the 152cm Loiano telescope 
(jCalchi Novati et al.ll2€07l ). This appears as a bright, AR ~ 19, short duration, ^fwhm ~ 3 days, flux variation, lying 
at about 3' from the M31 center, with a flat baseline on both previous OAB and POINT-AGAPE data sets. These 
characteristics, the excellent agreement with a Paczyhski shape together with the rather red color, i? — / ~ 1.1 (this 
feature is against the alternative hypothesis of contamination from some kind of eruptive variable, usually bluer), have 
all been con sidered as strong indications of the genuine microlensing origin of this flux variation. However, the analysis 
presented in ICalchi Novati et al.l (|2009af ) suffered from the serious limitation in that data points were available on the 
rising part of the flux variation only. This did not allow us to test the expected symmetric shape characteristic of 

^ http://archive.noao.edu. 

* The KPNO and OAB images have pixel scales of 0.26" and 0.58", respectively. 
^ P. Stetson, private communication. 

® We recall that the F814W filter is similar to the standard Cousin-/. 
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Fig. 1. — OAB-N2 light curve and residuals with respect to the best fit. Circles (boxes) are for R (/)-band, with filled (empty) symbols 
for OAB (WeCAPP) data. Top and middle panels; the solid (dashed) lines show the best fit with (without) inclusion of finite size source 
effect. The abscissa units are time in days (JD-2450000.0). The ordinate units are flux in ADU s^^ per superpixel. 



TABLE 1 

OAB-N2 : CHARACTERISTICS AND PaCZYNSKI FIT RESULTS (CORRECTED FOR FINITE SOURCE SIZE AND LIMB DARKENING EFFECTS) 





to (JD-2450000.0) 




p iE (days) 




R-I 


0''42'°50=.31 41°18'40".l 326.16 


4466.23loio 


0.07+0-07 


n 1^4+0-062 a 1+2.7 
U--'-34_o.o43 ^■-'--2.1 


1.245^013 22.25 ±0.51 


1.24 ±0.03 





Note. — It results x^/dof = 1.24. The reported error are for Ax^ = 1.0 (68.3% level). For the impact parameter uq there is no 
minimum value at the requested confidence level (as shown also in Fig.[2]|. For the (normalized) angular radius p we find also a secondary 
minimum at p ~ 0.08 (correspondingly uq ~ 0.11 and t» ~ 0.55). The source magnitude and color values are not corrected for extinction. 



microlensing events and to definitively conclude on the microlensing nature of this flux variation. We are now able to 
overcome this problem as data poi nts along the descent of the bump have been made available to us from the WeCAPP 
collaboration (iRiffeser et al.|[200ll ). These data have been collected at the 0.8 m Wendelstein telescope in two pass 
bands, R and /, similar to those used for OAB data. 



The analysis of the joint OAB- WeCAPP light curve, shown in Fig. [T] together with the results of the fit to be 
discussed below, is extremely interesting. WeCAPP data alone, missing the peak of the flux variation, do not allow to 
detect this flux variation. They are however essential to constrain the flux variation as they nicely cover the descent 
and contribute to the already excellent coverage along the rising wing of the fiux v ariation, extremely important to 
robustly constrain the microlensing parameters (|Baltz fc Silkll2000t iDominikI |2009| ). The excellent agreement with 
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a Paczynski shape and the achromaticity are now probed along all of the flux variation. We confirm the genuine 
microlensing nature of 0AB-N2. This conclusion is enhanced by the fact that we are considering two completely 
independent data sets analyzed, moreover, following two different schem es of phot ometry r eduction, the "superpixel" 
photometry with empirical correction of seeing va riations for OAB data l Ansari et al. 199 7tlCalchi Novati et al.|l2002l) 



and difference image analysis for WeCAPP data ( Alard fc Luptonl ll998': 'Gossl fc Riffeser 2002). For a recent analysis 
on difference imagi ng photometry specific for high surface brightness targets such as the M31 bulge region we refer to 
iKerins et all (|2010D . 

The excellent sampling along the bump allows us to move further in the analysis and to look for a b etter charac- 
teriza tion of this microlensing event. In particular we want to look for signatures of finite source effect (jWitt fc Mad 
Il994f) with the purpose to eventually constrain the lens proper motion ()Gouldiri994[ ). 

The relevant parameter for the finite source effect is p = where 0* is the angular radius of the source and 

6e the angular Einstein radius. We recall that finite source effects become releva nt whenever p ge ts larger than the 
impact parameter u. For the event amplification, following the scheme outlined in lYoo et al.l ()2004r ). we write 

Ait)^-Ap,,,{t)Biiz)Tx + Ak. (1) 

Here Apacz is the standard Paczynski amplification, is the amplification including the finite source effect, for which 
we us e the exact iWitt fc Mad ()1994D expression (for a different approach we refer to the recent analysis in iLee et ahl 
|2009| ). z = u/p, a nd F^ are the, wa velength dependent, linear limb darkening coefficients. The expression for Bi{z) is 
given in Eq. 16 of lYoo et all (|2004[ ). 

The inclusion of limb darkening is relevant as we know the source flux n ot to be uniformly distributed across the 
star surface. To evaluate the Fa we make use of the results of Ivan Hammd ()1993l ). Starting from the event color (the 
driving and essential parameter) we have to further s pecify the star temperature, surface gravity and metalicity To 
this purpose we make use of the iMarigo et al.l (|2008l ) isochrones. We consider two cases: a 12 Gyr populatio n with 
Z = 0.030, typical for the bulge, and a 2 Gyr population with Z = Zq = 0.019, typical for the disk (Riffeser et al.l 
12006*). From the light curve analysis, and the magnitude zero point, we es timate a deredde ned color {R — I)c = 1.24. 
Taking into account the MW foreground extinction, E{B-V) = 0.062 ( Schlegel et al.lll998[) . we get {R~I)a,c = 1-20. 
For bulge/disc sources we then get Tes ~ 3600/3550 K, log(5) - 0.90/0.60. Finally we get F = 0.81/0.81 (6.62/0.60) 
for Rc [Ic], respectively. Even if we include a term of internal extinction (jRiffeser et al.l 120061 ) the final result does 
not change significantly. Using ext^ = 0.19 for the bulge we finally get to F = 0.80 (0.59), and with E{B — V)= 0.22 
for the disc we get to evaluate F = 0.76 (0.58 ). In conclusi on, we use F^ = 0.80, 0.60 for R and / band data. We reach 
consistent results basing our analysis on the IClaretl ()2000[ ) limb darkening coefficients instead. 

We carry out the light curve Paczynski fit (corrected for finite source size and limb darkening effects) using Eq. [T]for 
the microlensing amplification*. For a simultaneous fit of the two joint data sets with two bands each, overall we have 
12 independent parameters (4 and 8 to account for geometry and fluxes, respectively). The microlensing amplification 
parameters: the Einstein time, ^e, the minimum value of the impact parameter, ug, the (normalized) source angular 
radius, p and the time at maximum amplification, ig- The 8 fiux parameters are, for both R and / and OAB and 
WeCAPP data, the background value and the source flux, (fiR and 0/. 

Our first purpose is to constrain the normalized source angular radius, p, or equivalently = pt^, the source radius 
crossing time (the Einstein time turns out to be better constrained with respect to p). Then, since the event source 
fiux and color can be used to estimate the source angular radius, 6^, to finally constrain /i, the lens proper motion 

p, = — = — (2) 
tB 

As a tool we make use of a analysis in which we scan the parameter space of one (or more) fit parameters while 
the remaining are left free and determined, at each step of the scan, by the fit^. The main results of this x^-based 
analysis, discussed below, are presented in Table [1] and Fig. I2I3I 

A first relevant outcome is the excellent agreement with a Paczynski shape refiected in t he small value / dof = 1.24 . 
The best fit value for the source magnitude is significantly fainter than in our previous iCalchi Novati et al.l ()2009aD 
an alysis, although affec ted by a very large error. In particular, the result is 1.3 mag fainter (beyond 95% level) than 
the lMassev et al.l ()2006f ) catalog source. Furthermore, it is in agreement with the HST-based estimate in Sect. 12. ll (for 
/-band magnitude 21.0 ± 0.5 to be compared with 20.8 ± 0.3). Both the source magnitude and the Einstein time are 
rather well constrained within the fit: at 95% level we obtain i?* = 22.215'^ and t^ = 8.2I3 g, respectively. This is not 
the case, however, for the normalized angular radius p, which is degenerate with the impact parameter uq (the impact 
parameter is in fact unbounded from below even at the 68.3% level) along two different directions (Fig. [5J top panel). 
Speciflcally, we find an absolute minimum at p w 0.15 with uq ~ 0.07 (along the corresponding direction of degeneracy, 
Uq goes to zero with p > 0.1) and then a secondary minimum, with A^^ ~ 0.9, with p w 0.08 and uq « 0.11 (here, just 
above the 68.3% threshold for Ax^ = 1.0, p goes to zero with roughly constant uq). The former rather large value of 

We make use of the relationship F = 2u/ (3 — u ) to lin k F to the linear limb darkening parameter u (not to be confused with the 
microlensing impact parameter) given in Ivan Hammel II1993I) . 

* We rescale the error bars, by ~ 10%, so to force, along the baseline, x^/dof to unity. This only marginally, however, changes the 
outcome of the analysis. 

^ We use, within the CERNLIB MINUIT package, the minimizatio n tool MIGRAD "a variable-metric method with inexact line searc h, 
a stable metric updating scheme, and checks for positive-definiteness" |http : / / wwwasdoc . web . cern . ch / wwwasdoc /minuit/minmain.html[ 
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Fig. 2. — 68.3%, 90% and 95% contour maps in the parameter space p-«o (top panel) and p-t^,- p, uq and are the normaUzed 
angular radius of the source, the impact parameter and the Einstein time, respectively. 

p, with p/uq K,2 might be taken as an indication of finite size source effect. Given the caveats further addressed below 
this would make of 0AB-N2 the first microlensing event reported toward M31 with a signature of finite size source 
effect. Indeed, the secondary minimum, as well as the full degenerate parameter space region for p < 0.1, just above 
the Ax^ = 1.0 threshold, do not allow us to draw firm conclusions on this issue. Overall, the indication for p > 
remains extremely tenuous. In particular the improvement with respect to the case with no inclusion of the finite 
size parameter, p, is only marginal, Ax^ ~ 1. The inclusion of p, however, is essential for the analysis as a tool to 
allows us to constrain the lens proper motion. The degeneracy for p is shown also in the parameter space p-t^ (Fig- [3 
bottom panel). This result is refiected in the estimate of the source radius crossing time for which we find the absolute 
minimum at « 1.25 days, but also a secondary minimum at t^, « 0.55 days (Fig. |21 top panel). 

The missing information to get to the lens proper motion p is the s ource angular ra dius, 0,. To this purpose we use 
the empirical relationship ^ = 9if(V.V — K ) for giants discussed in van Bellel (|1999D . Taking into account the MW 
foreground extinction (Schl egel et al.llTQQSl ). and the Marigo ctaL, (,20081 ) isochrones as above we estimate, both for a 
bulge or a disc source, V — 23.3 and V — K = 5.0. This gives 

0, -0.67t°:il pas, (3) 

where the error budget is dominated by the statistical uncertainty of the source magnitude within the light curve fit. 
The result we get to for 6'* is very similar, within 1%, if we additionally take into account the internal M31 bulge or 
disc extinction as discussed above. In these cases we would find a fainter source with a bluer color, and these two 
effects balance one each other in the evaluation of the angular radiu s. (An alternative, more recent and more accurate 
0* — 0*(V^, V — K) relationship is discussed in lKervella et al.ll2004 The result we would get to in th is case for 0^, is 
about 20% larger, still comfortably within la. For our discussion we however prefer that discussed in Ivan Belle T99l 
as it is derived from a sample of data with 2.0 <V — K < 8.0 whereas iKervella et al.|[2004l use a more Cepheid-specific 
sample with 1.0 <V~K< 2.4, with the former being more suitable for us given the source color). For the discussion 
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Fig. 3. — Plots of vs the source radius crossing ti me (t t, top panel) and vs the lens proper motion (/x). Bottom panel: we use 

The solid, dashed and thin dashed lines indicate the 68.3%, 90% and 95% 



an estimate of the source angular radius of 0.67 fj-as (Sect. I2T21 
confidence level, respectively. 

of our results on /i we finally recall the useful relationship 



(4) 



where the coefficient a depends only on the color (0 is the source flux). In particular this allows us to carry out the 
analysis (Fig. |31 bottom panel), in the 3-d parameter space ^e, P and (/). 
The value for the lens proper motion we obtain, corresponding to the best fit value, is fi — 0.93kms-ikpc~\ The 
discussed large residual degeneracy for the angular radius p, and the related aspects discussed above, lead us however 
to report, as a more robust outcome of the present analysis, an upper limit for p (and t^,) and a lower limit for fj, only. 
At 68.3% level we get 

p<0.22, 

< 1.35 days, 
^>0.75 kms^^kpc"^. 

For the lens proper motion we also get fi > 0.68 (0.64) kms^^ kpc^"^ at the 90% (95%) level, respectively. 

3. DISCUSSION 



The results of the light curve analysis presented in Sect. 12.21 that in particular confirm the outcome on the source 
magnitude discussed in Sect. 12. 1[ can be used to establish the lens nature of the 0AB-N2 event. As previously 
discussed, the relevant issue is to distinguish between self lensing and MACHO lensing. To this end we focus on the 
lens-source relative proper motion, fi, because it is independent of the lens mass and so depends only on the kinematics 
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Fig. 4. — Expected distributions for the lens proper motion, fj,, along the line of sight toward the OAB N2 position for different source-lens 
populations evaluated within the Monte Carlo simulation of the experiment. Top panel: thin solid, dashed and dotted lines are for bulge 
sources and bulge and disc lenses and disc-disc events, respectively; thick solid line for the overall M31 self lensing (each distribution is 
normalized to 1 separately). Bulge- bulge events are expected to contribute for about 80% of the expected signal; MW disc events (not 
shown, for which the lens proper motion distribution is expected to peak well beyond 1 kms"'^ kpc"'^) are expected to be negligibile (see 
text for details). Middle (bottom) panels: M31 (MW) MACHO lenses with the overall (M31-I-MW) MACHO distributions normaUzed to 
1, and 75% signal expected from M31 lenses (in the inset, bottom panel, showing a zoom for the MW MACHO distribution, the units 
are 10"'^). The solid, dashed and thin dashed vertical lines (68.3%, 90% and 95% level, respectively) indicate the lower limit for the lens 
proper motion derived in the light curve analysis (Sect. [22] Fig. [Sj. 

and spatial distribution of the lens population. The po t ential role played by t he len s proper motion, also in relation 
to M31 pixel lensing, has been first discussed in iGouldl (|1994[ ): iHan fc Gouldl ()1996[) . In these analyses the emphasis 
was given in particular to the possibility to constrain MW MACHO lenses with respect to M31 lenses. This is relevant 
as MW MACHOs are expected to contribute, for given MACHO mass and mass halo fraction, for about one third of 
the overall expected MACHO signal. The identification of MW lenses would therefore give a strong indication for a 
MACHO signal. In the following we will point out how and to which extent an analysis of the lens proper motion may 
allow one to distinguish also between the different M31 lens populations. 

As a tool of analysis, to compare wit h the outcome of the light c urve fit, wc will make use of the Monte Carlo 
simulation of the experiment described in ICalchi Novati et al.l ()2009al) . As a prior we consider the line of sight toward 
the 0AB-N2 position (this gives however only marginal quantitative differences from the more general case in which 
the overall PLAN field of view is considered). We also recall that as an output we have only light curves that are 
selected within the Monte Carlo (where in particular the experimental set up is carefully reproduced together with a 
basic selection pipeline). 

At the 0AB-N2 position, only « 3' from the M31 center, the dominant contribution from M31 self lensing is 
expected from events with bulge sources and bulge lenses. Specifically, we expect ^ 80% of events to belong to these 
populations, with the rest equally divided between bulge-disc (both source-lens and lens-source) events, with only a 
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marginal 1% of disc-disc events. In Fig. |4]we present the results of the Monte Carlo simulation for the lens proper 
motion distributions for different source and lens populations. The most striking outcome for our purposes is the 
threshold at about /ith ^ 1 kms~^ kpc""'^, such that for /i > /ith we expect almost no M31 lens events, while on the 
other hand almost all of MW lens events are expected to fall in this region. This is easily explained on the basis that 
finite size effect is never important for MW events (for which on average the Einstein angular radius, proportional 
to the square root of the lens-source distance over the lens distance, is about one hundred times larger than for M31 
events and p oc 9^^) and for which, on the other hand, the lens distance is by far smaller. The relatively large lower 
limit value we find for the lens proper motion might therefore be taken as an indication for the lens to belong to the 
MW halo. (This outcome, however, would be clearly at odd with the discussed marginal indication we find for finite 
size source effect.) However large, the lower limit for the lens proper motion is still fully compatible with M31 lensing. 
In particular, among M31 lensing events, we find that a broader distribution for /i is expected for MACHO lensing, 
while the limit we find is hardly compatible with the expected distribution for bulge-bulge events. Although we can 
not certainly rule out, for this specific event, the self-lensing hypothesis, this outcome strongly suggests the MACHO 
lensing hypothesis for the 0AB-N2 event. 

More specifically, corresponding to the 68.3% lower threshold obtained above, /i > 0.75 kms""'^ kpc~^, for M31 
MACHO lensing we still have about (for bulge source events and similarly for disc source events that represent, 
however, only about 15% of the expected M31 MACHO signal) 6% of the expected signal (10% and 14% respectively 
at the 90% and 95% levels). On the other hand, at the 68.3% level, we expect less than 1% of the overall bulge-bulge 
signal (still about 70% of the overall M31 self lensing signal) and about 1% of bulge-disc and disc-bulge lenses. At 
95% level these ratio rise to almost 2% for bulge-bulge events (about 60% of the overall M31 self lensing signal) and 
to almost 5% for bulge-disc and disc-bulge lenses. Although in an absolute sense M31 disc-disc events give only a 
marginal contribution, it is worth noticing that, along this M31 central line of sight, they show a marked bimodal 
distribution for the lens proper motion with a long tail for large values of ^. In particular, larger than the 68.3%, 
90% and 95% threshold values we still find, respectively, 2.6%, 7.2% and 12% of the expected events belonging to this 
population. Overall, at the 68.3%, 90% and 95% levels, we expect respectively about 0.6%, 1.3% and 2.1% of the M31 
self-lensing signal. 

Besides M31 bu lge and disc lenses, in principle one might cons ider also Galactic disc lenses. Using standard models 
for the MW disc (|Han fc Gouldl[2003L ICalchi Novati et all 120081 ) the ratio of the expected number of events for this 
population is below 1% with respect to M31 self lensing^°. This res ult, joined to the overall ou tcome of one, possibly 
two, microlensing events reported from the 2007 PLAN campaign (jCalchi Novati et al.ll2009"al) makes the possibility 
for 0AB-N2 to be a MW disc event extremely unlikely. However, for MW lensing, the lens proper motion is much 
larger than for M31 lensing and this applies even more strongly to stellar lensing than MACHO lensing for these two 
galaxies. In fact, typical values are well above 1 kms^^kpc"^, larger, on average, than for MW MACHOs. Hence, 
in particular almost all of the expected MW disc events fulfill the lower limit condition we find on the lens proper 
motion. Larger than the 95% level lower limit, the even small fraction of expected M31 self lensing events is however 
still larger than the overall expected MW disc signal. The overall small number of expected MW disc events makes 
therefore this contribution always negligible (so that in particular we do not show it in Fig. 14]). 

As for the Einstein time scale, which is rather well constrained from the light curve fit, comparing to the Monte 
Carlo distributions we find the fit result to match the expected self-lensing signal characteristics well. For MACHO 
lensing the result is dependent on the lens mass. The fit result is fully compatible with MACHO in the mass range 
(10~2 - 10"^) M©, but still compatible with even larger MACHO mass as 0.5 M©. 

Finally, we comment on our result for the finite source size effect, with the best fit value indicating a signature for 
p > 0. As already pointed out this effect is expected for M31 lenses only. Still, the value we find (as reflected also in 
the proper motion best fit value) has an a priori low probability. Together with the only marginal detection, with a 
Ax^ ~ 1 to separate from the case with no finite size effect, this makes this detection only "possible". 

4. CONCLUSION 

In the present paper we have presented an updated a nalysis o f the M 31 microlensing candidate event 0AB-N2 
first reported by the PLAN collaboration in Calchi Nova tTet al.l ([20093) ■ In this previous work, in particular, the 
insufficient sampling along the flux variation did not allow us to probe unambiguously the microlensing nature of this 
flux variation. Here we have taken advantage of newly available WeCAPP data to complete the light curve sampling. 
From a methodological perspective, we observe that the importance of merging different data sets for microlensing 
observation might look even trivial when considering that Galactic Bulge events are routinely observed by multiple 
telescopes when searching for planetary anomalies. The present analysis, however, shows the extent to which this is 
suitable, feasible and, as this case nicely shows, possibly essential also for M31 pixel lensing events. (Such a strategy, 
besides, is unavoida.ble in order to achieve the challen ging pu rpose of a systematic research for exoplanets in M31, 
ICovone et alll2000t iBaltz fc Gondolol 120011: IChung et a l. 2006 : llngrosso et all [20091 A first attempt in this se nse is 
being carried out by the ANGRSTOM collaboration [Kerins et al l 120061: iDarnlev et all 120071: iKim et al.l [20071 The 
usually short expected timescale, ipwHM, makes however potentially important such a strategy also for "normal" 
single-lens events.). 

A preliminary analysis on archival images of the 4m Mayall-KPNO telescope and of WFPC2/HST allowed us to 

Furthermore, a sizeable fraction of these MW events is expected, because of the extremely small observer-lens distance, to have too 
bright lenses for being suitable candidates. 
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reject the tentative identification of the 0AB-N2 source made in ICalchi Novati et al.l ()2009al ) with a rather bright 
iMassev et alJ (|2006l ) M31 catalog source. The event source, in fact, appears to be a nearby fainter blended object on 
the KPNO image, as confirmed also by the HST analysis. Although we can not conclusively determine the source 
magnitude, because of blending and the huge background level, even this partial result is a first step toward a better 
characterization of this event. Indeed, the new light curve photometry analysis confirms that the source must be a 
fainter object than previously concluded. 

The improved sampling along the 0AB-N2 light curve then allows us to reach a few important conclusions. First, 
we find a strong confirmation of the microlensing nature of this fiux variation. The fiux variation characteristics, the 
large signal to noise ratio, the short duration and the position, within the inner 3' of M31 where the expected rate 
is larger, together with the further evidence provided by the joint analysis of two completely independent data sets, 
put this fiux variation among the most convincing and interesting M31 microlensing events reported so far. Second, 
we have carried out a more detailed analysis in the event microlensing parameter space with the specific purpose to 
constrain the lens nature, whether due to self lensing or to MACHO lensing. The main outcome we get to is a lower 
limit for the relative source-lens proper motion, /i. As we have shown, this parameter, independent of the lens mass 
and therefore less subject to model issues, not only may allow one to distinguish M31 lensing from MW lensing (where 
the main expected signal is MACHO lensing and for which /i is significantly larger than for M31 lensing) but also 
among different M31 lens populations. Specifically, at 95% level, we find fi > 0.64 kms^^kpc^ . Comparing to the 
output of the Monte Carlo simulation of the experiment we conclude that, if not to the MW halo population, this 
value is more compatible with M31 MACHO lensing rather than with self lensing, for which only about 2% of the 
overall population is expected to fulfill this constraint. Although the self lensing hypothesis can not be definitively 
ruled out for this single specific event so that we are not yet in a position to strongly discriminate self lensing from 
MACHO lensing, this analysis shows the extent to which the lens proper motion, with no other prior assumptions on 
MACHOs but their spatial distribution and kinematics, is a powerful tool toward this purpose. 

As an intermediate step of the light curve fit x^-based analysis we have looked for possible signatures of finite 
size source effect (through a study of the normalized angular radius, p) for which however we find only marginal 
indications. We recall that an evidence in this sense would imply that the lens belongs to some M31 populations. The 
robust outcome we get to is rather an upper limit for p and on the related quantity, the source radius crossing time, 
t* = pt^ which, together with an estimate, from the source magnitude and color, of the physical angular radius, finally 
gives us the aforementioned lower limit for the lens proper motion. 

The outcome of this analysis is affected by the lack of an accurate and precise estimate of the unresolved source 
magnitude and color, independent of the light curve fit, for which a focused proposal on HST data would be extremely 
welcomed. 

In the frame work of M31 pixel l ensing, the present anal ysis of 0AB-N 2 is the third example, following those for the 
events PA-Nl (jAuriere et al.ll2001[ ) and PA-S3/GL1 (Riffe ser et al.ll2008l ). where a detailed investigation of the lensing 
parameter space allowed to draw conclusions on the lens na ture. Inte restingly, although in the present work with 
a worse confidence level with respect to that reported in iRiffeser et al.l (2008 .) for PA-S3/GL1 and with the caveats 
discussed in the Introduction as for PA-Nl, in all cases the analyses point toward the same outcome, namely that 
MACHO lensing should be preferred over self lensing. 

We warmly acknowledge the WeCAPP collaboration, in particular A. Riffeser and S. Seitz, for making available 
to us, prior to their pubhcation, the WeCAPP data relative to the 0AB-N2 light curve. SCN, VB, LM and GS 
acknowledge support by MIUR through PRIN Prot. 2008NR3EBK_002, and by FARB of the Salerno University. 
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